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Abstract. In this paper we investigate spectral properties of Lapla- 
cians on Rooms and Passages domains. In the first part, we use Dirichlet- 
Neumann bracketing techniques to show that for the Neumann Lapla- 
cian in certain Rooms and Passages domains the second term of the 
asymptotic expansion of the counting function is of order For the 
Dirichlet Laplacian our methods only give an upper estimate of the 
form vA- In the second part of the paper, we consider the relation- 
ship between Neumann Laplacians on Rooms and Passages domains 
and Sturm-Liouville operators on the skeleton. 



1. Introduction 

Let — Ajv ; n denote the Neumann Laplacian on a bounded open subset 
of M. n ,n > 2, and let E : H l {VL) — > L 2 (Vt) be the canonical embedding, 
where H l {yi) is the standard Sobolov space. Then (— Ajv,o + = EE* , 
and — Ajv,n has a compact resolvent (and thus a discrete spectrum) if and 
only if E is compact. The compactness, or otherwise, of E is determined 
by the nature of the boundary dVt of f2. In pQ it was shown that E be- 
ing compact is equivalent to a generalised extension property, which is that 
there exists a function space A(M. n ) on W 1 which is compactly embedded in 
L 2 (B) for any ball B C W 1 , and is such that there is a continuous extension 
£ : H 1 ^) — > A(W n ). In general A(R n ) is a space of arbitrary smoothness. If 
<9f2 E Lip a ,0 < a < 1, i.e., dCl coincides with the subgraph of a Lip a func- 
tion in a neighbourhood of each point, then A(M n ) = H a (M n ), the Sobolev 
space of order a, so that, in this case, there is a reduction of smoothness in 
going from Q to M. n . For domains with singular boundaries, like the "Rooms 
and Passages" domain defined in Section 2, the compactness of the embed- 
ding E can be taken as a measure of the smoothness of the boundary, which 
is otherwise difficult to describe. When E is not compact, [7] shows that for 
any closed subset S of the non-negative real numbers, there exists a mod- 
ified "Rooms and Passages" or a "Comb" domain such that S equals the 
essential spectrum of —Anq- 



The authors would like to thank the referee for carefully reading an earlier version of 
the manuscript and Dr. Rob Davies for producing the figures. 

1 



2 



B. M. BROWN, W. D. EVANS, AND I. G. WOOD 



Suppose that — Ajv.n has a discrete spectrum and denote the number of 
its eigenvalues less than A by Nn(X). The problem of determining the as- 
ymptotic behaviour of N]y(X) as A — > oo has a long history. If the boundary 
d£l of fi is sufficiently smooth, it has long been known that one has the 
asymptotic formula 

(1.1) N N (X) = cj n (2ir)- n \n\\ n / 2 + o(A n / 2 ), 

where uo n is the volume of the unit ball in W 1 and \£l\ is the volume of the 
domain Q. For general f2 the Weyl term uj n {2-K)^ n \Q,\X n ^ 2 dominates, and in 
the error R(X) := N N (X) - w„(27r)- n |fi|A n / 2 , it is the boundary dtt, rather 
than any other topological or geometrical feature of which is dominant. 
For instance, in [9], it is shown that if d£l G Lip a ,0 < a < 1, the remainder 
R(X) = 0(X^ n ~ 1 ^ 2a ) and this is order sharp. Here and in the following, 
the 0(/(A))- and o(/(A))-notation is to be understood in the limit A — > oo. 
Moreover, in [3], a domain of von-Koch snowflake type is exhibited for 
which R(X) x X d °, where d Q denotes the outer Minkowski dimension of the 
boundary and ^4(A) x -B(A) means that |^4(A)/i?(A)| is bounded above and 
below by positive constants as A — >• oo. Specific information about the error 
is hard to come by. For general domains one of the few available techniques is 
the Courant-Weyl variational method involving "Dirichlet-Neumann brack- 
eting" . This effectively reduces the problem to estimating the corresponding 
counting functions Nd{X) and Nn(X) for Dirichlet and Neumann Laplacians 
on a set of cubes which cover f2. A variant of this technique is used in [9] 
in which cubes are replaced by other relatively simple sets for which the 
Dirichlet and Neumann Laplacian counting functions can be estimated from 
above and below. The first part of this paper is a contribution to the study 
of how the error term R(X) depends on the boundary d£l. We look in detail 
at the much studied "Rooms and Passages" domain, in which the Weyl for- 
mula (II. ip holds, being particularly concerned with upper and lower bounds 
for R(X)X- 1 / 2 . 

In [5] it was shown that for a rather restricted class of domains f2, (which 
does not include Rooms and Passages) — Ajv,n has a compact resolvent if and 
only if a Sturm-Liouville operator defined on the skeleton of has a compact 
resolvent. Recall that the skeleton of an open set f2 is the complement of 
the set of points x in Q for which there exists a unique point y on (90 such 
that |y — x| is equal to the distance of x from M n \ fi. The result in [5] was 
motivated by Theorem 3.3 in [2] in which is a horn, whose skeleton is a 
half-line. In the second part of the paper (Sections 4 and 5) we investigate 
this problem for general Rooms and Passages domains. 

2. Rooms & Passages domains 

We consider a Rooms and Passages (R&P for short) domain Q defined as 
the union of square rooms Ri {i odd) of size hi x hi joined by rectangular 
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passages Pj (i even) of size hi x <5j with 5, < min{/ij_i, /ij+i}. We assume 
hi — >• as i — > oo. See Figured! 
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Figure 1. The start of a general Rooms and Passages domain. 

In this section and the next, we further restrict the R&P domains under 
consideration by assuming that hi = C l and 5i = kC ta , where C < 1, a > 1 
and k are constants, with k < C 3 ~ 2q . This guarantees that each passage is 
narrower than both adjacent rooms. It follows from the general result in |3j, 
Example 6.1.1, that the embedding E : H l {Q) — > L 2 (Q) is compact if and 
only if a < 3. 

Proposition 2.1. For a > 3 we have 6 o~ ess (— A^q) where a ess denotes 
the essential spectrum of the operator. 

Proof. We give a simple singular sequence proof. By Qj we denote the 
domain consisting of the first j rooms and passages while Tj = ft \ Qj 
denotes the "tail". Choose a sequence of C°° functions ipj which depend 
only on the x- variable such that supp cpj C Q±j \ 



(2.1) 



on 
on 



n 



4j-l 



\n 



2j, 



% u n 2j . 



and 



o(c-^) on n 2j \n 2j -i, 
-i, - ^ OiC on n 4j \n 4 j_i. 

In other words, the function <pj is initially zero, increases to 1 in the j- 
th passage and falls back to zero in the 2j-th passage. The support of 
the sequence of the <pj 'disappears' into the tail and each (pj satisfies the 
Neumann boundary condition on d£l. The singular sequence we wish to 
consider is given by 

(2.2) fj(x,y) = T . — 7T <pj(x) ior(x,y)en. 



Obviously, \\fj\\ L 2 {n) 
We have the following: 

( 2 - 3 ) ll^'llx 2 (n) 



1, while fj weakly converges to in L 2 (fi) as j 



oo. 



i-.ii 



1/2 



as j 



oo 
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and 

4j — l 4j 
1=1, j odd i=l, i even 

Therefore, 

|Oy|-|^| = C 2 — - 4 +kC 2 ^ _ c2{l+a) x as j — )• oo, 

and for large j, \\<Pj\\ L 2^ behaves like C 2j . 

As ip'j is supported only on two passages we get from (|2.ip . 
(2.5) 

IK'fe(fi) = 0(C- 8j kC 4j{1+a) +C~ 4j kC 2j{1+a) ) = 0{C~ 2j{l ~ a) ) as j -> oo. 

Hence, as j — » oo, 
(2.6) 

if a > 3. This implies that (— A) 1 ' 2 /.,- — >• as j — >• oo although H/jH^cj) = 1- 

Therefore, we have G cr ess ((— Aatq) 1//2 ), which implies E <7 ess (— A^n) 
when a > 3. □ 

3. Eigenvalue asymptotics of the Dirichlet and Neumann 

Laplacians 

The special class of R&P domains introduced in the previous section will 
be considered. Let a < 3 and denote by Njj(X), iV/v(A), respectively, the 
counting functions of the Dirichlet and Neumann Laplacians on the R&P 
domain f2. To determine their asymptotic behaviour, we shall apply the 
Dirichlet-Neumann bracketing technique. 

We begin this analysis of the spectrum of the Laplacian on an R&P 
domain 0, by discussing the contribution from the tail. As before, let 
Q = 0,2m U T2M where f^M consists of the first 2M rooms and passages 
and T2M denotes the tail. On applying Theorems 4.6 and 5.1 to Example 
6.1.1 in [3], it follows^ that the optimal constant K (T2M) in the Poincare 
inequality 

11/ " H m \\l*(t 2M ) <K(T 2M ) ||V/|| l2(3W) , (/ e H\T 2M )), 
satisfies 

K (T 2M ) < cC^ M 



^By Theorem 4.6, the bound on K(T2m) involves a quantity c(Ji), which according to 
Theorem 5.1 is bounded by Noting that their constant C corresponds to our C~ , 

Theorem 6.1 and Example 6.1.1 with k = 2M then give N a x C 2M(3_a) . 
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where c is a positive constant. 

It follows that the Neumann Laplacian on T 2 m will only contribute the 
trivial eigenvalue if 

p.!, »«i/c)»c«^, u. M> 2(3 ': g XsU 

Since a > 1, the tail after M passages has area 

I T 2M ]= ^ + kC \_ c g 1+a) * C" as M -> co. 

It follows that for M as in (I3TT]) 

|T 2M |=0(A- 2 /( 3 - a )) 
and so _^ 

= |JW| +o(v / a). 

Now, for M satisfying (|3.ip . the counting function for the Neumann- 
Laplacian on differs by at most 1 from that on £l2M- It follows that 

(3.2) R(\) - R 2M (X) = o(\/A), 

where R(X) and R 2 m{^) denote the error term for the problem on O and 
^2A/, respectively. 

As we have that the (re + l)-th Neumann eigenvalue is a lower bound 
for the re-th Dirichlet eigenvalue (see [6]), for A satisfying (|3.ip . the tail 
cannot contribute any Dirichlet eigenvalue, so the same reasoning as for the 
Neumann case implies that (13. 2p also holds in the Dirichlet case. 

3.1. Asymptotics for iVjv(A). Our strategy here is to partition the domain 
and use the Dirichlet Neumann bracketing technique to obtain the required 
estimates. In order to obtain these estimates we first obtain a lower bound 
for the number of eigenvalues of the Neumann-Laplacian, we partition the 
rooms into five sections imposing Neumann boundary conditions on the 
boundary of O and Dirichlet boundary conditions on all artificially intro- 
duced internal boundaries (see Figure [2]). An upper bound is obtained by 
only introducing an artificial boundary to separate the room from the neigh- 
bouring passages and imposing Neumann conditions on all the boundaries 
(see Figure [3]). This is a simple consequence of the variational principle. 
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Figure 2. Artificial boundaries and sub-regions of a room 
for the lower bound: Neumann conditions on the dotted lines. 
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Figure 3. Artificial boundaries of a room for the upper 
bound: Dirichlet conditions on the dotted lines. 



We first consider the one-dimensional problem on [—a, a] and have the 
following eigenfunctions and eigenvalues: 

(a) Dirichlet conditions at both end points: 

rmr(x + a) 



csm ■ 



A; 



m 2 7r 2 



2a 4a 2 ' 

(b) Dirichlet conditions at —a and Neumann at a: 



m > 1. 



csm • 



(2m + 1)tt(x + a) 



A, 



(2m + l) 2 ^ 2 



4a 7 16a 2 

(c) Neumann conditions at both end points: 



m > 0. 



m7r(a; + a) 
^ m (a;) = ccos - , A m 



m 2 7r 2 



m > 0. 



2a ' 4a 2 
For the lower estimate for a room, this leads to the following. 
(1) In regions I and V, we have the set of eigenvalues 



m 2 vr 2 (2n + 1) 2 tt 2 
A m , n = , o H : m, n > 



with a = — , 6 

2 ' 



4a 2 166 2 

(2) In regions II and IV, we have the set of eigenvalues 

f> _ (2m+l) 2 7r 2 n 2^2 "I 
1 Am,n — Jg^2 r -^r ■ m d. u, n ^ L > 

with a = 6 = | (C ^- 1 ) - C a( - j+1 ^) . 

(3) In region III, we have the set of eigenvalues 



Am,n 



m 2 7r 2 ra 2 7r 2 



+ 



4a 2 46 2 



: m, n > 1 



C«Cf+i) 

with a = — , 6 = /c 

2 ' 2 



We now need to count the integer lattice points in the first quadrant 
satisfying \ m , n < A. By [8], the number of integer lattice points in a plane 
region X ■ R, where X is a real scaling parameter, is given by 

. 131 



Af(X) = AX 2 + o(Xi 



as X — y oo 



for any e > 0, where A is the area of the region R. This gives us the 
following: 
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(1) In regions I and V, we have an ellipse with half axes 2a ^ and 2b ^ 
which is shifted by —1/2 in the y-direction. As an upper estimate 
for the area below the x-axis we use the area of the rectangle which 
is subtracted from the area of the quarter ellipse (the error from 
this can easily be seen to be o(\/A)). We then add the contributions 
along the two axes to get 

, f\ m 2 7r 2 (2n + l) 2 ir 2 
card <^ X m , n = 4q2 H : rn,n>0, A m , n < A 

ir2a>/\2b>/\ 12a\f\ . r-, 2aV\ 2bV\ 
> o vA H 1 

4 7T 7T 2 7T 7T 7T 

= ab x+ (2b + a)VX_ ^ 

7T 7T 

Here, card A denotes the cardinality of the finite set A. 

(2) In regions II and IV, we again have an ellipse with half axes 2a ^ 

and 2b ^ which this time is shifted by —1/2 in the x-direction. As 
an upper estimate for the area left of the y-axis we use the area of 
the rectangle which is subtracted from the total area. We then add 
the contributions along the y-axis to get 

f (2m + l) 2 7r 2 n 2 vr 2 
card <^ \ rn . n = h : rn>0, n>l, A m , n < A 

ir2ay/\2by/\ \2bJ\ , - 26v/A 
> o VA H 

4 7T 7T 2 7T 7T 

= ^A + ^-o(VA). 

7T 7T 

(3) In region III, we count the lattice points in an ellipse with half axes 
^ and to give 

card j A m , n = + : rn,n>l, A m>n < A 

7T 2avA 2&vA , a& , , /-, 

> -— o VA = — A — o vA . 

4 7T 7T 7T 

Collecting these results, we obtain a lower estimate for the contribution of 
the j-th room to the counting function of the form 



C 2 i 

Nj > —A + 

J 47T 



k / 



^ _ 0( VA). 

7T 



For the upper estimate of the counting function, we simply need to con- 
sider the eigenvalues 



m 2 7r 2 n 2 7r 2 
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of the Neumann-Laplacian on the square, where a = b = Of. In order to 
count the integer lattice points, we take the area of the ellipse with half axes 
2a ^ and 2b ^ and add the additional points along the x- and y-axes: 

f m 2 ir 2 n 2 7r 2 

card < A m , n = ^ 2 + : m, n > 0, A mj „ < A 

vr2a\/A2&\/A 2a^/\ 2b\f\ , ^ 
< + + + o(VA 

4 7T 7T 7T 7T 

ab x 2aV\ 2bVX , r-. 
= — A + + + o(VA). 

7T 7T 7T 

As an upper estimate for the contribution of the j'-th room, we therefore get 

(~<2j r)f<j 

Nj < ^— A + — Vx + o(V\). 

J 47T 7T 

The calculations involving the first room are a little different and we only 
use three partitions (see Figure 2]) . 




Figure 4. Subdivisions of the first room for the lower 
bound: Neumann conditions on the dotted lines. 



'1) In regions I and III, we have the set of eigenvalues 

m 2 TT 2 (2ra + l)V I , C , C-kC 2a 
1^ + 166 2 : m ' n "° With a= 2' b = 4 ' 



(2) In region II, we have the set of eigenvalues 

J (2m + 1) V ^ n 2 TT 2 
\ 16a 2 



+ 



q kC 2a 

,,,, : rn > 0, n > 1 }> with a = — , b = — - — . 
4o z 2 2 



Proceeding as for the other rooms we again need to count integer lattice 
points in the first quadrant. 

• In regions I and III, we have 

. f\ m 2 TT 2 (2ra+l) 2 vr 2 
card <( A m , n = - 9 H ^2 : n > 0, A m , n < A 

o(V\). 



4a 2 

aft, (2b + a)V\ 



> — A + 

7T 



7T 
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• In region II, we have 

f (2m + l) 2 7r 2 n 2 7r 2 
card < A m , n = h -^p- : m > 0, n > 1, A m , n < A 

^ a& , 6a/A , 

> — A + o(VA). 

7T 7T 

Combining these we get a lower estimate for the contribution of the first 
room 

We next sum over the rooms (omitting the o(\/A)-term for convenience). 
Denoting the volume of the first M rooms by V(Rm), we see 

2M 2M AC , 2 j 



£ ^ a £ 



47T 

j=l,j odd j=l,j odd 



+ [2 ^ C i-^ic 2a + ]T ^(l + C 2 ") ^ 

\ i=l J odd \ i=3,i odd / / 

" 4vr A + ^ 2C 1 - C 2 2 ' 1 - C 2 <* J vr ' 

Also, 

^ AT ^ /AC 2 -? 2VXcA F(i? M ), 2V\ ( C 2 - C 2 ^)' 

Z-^i 3 — I + - I - ~^~ X+ ^C I ^2 

j=l,j odd j=i,j odd \ / \ 

We now proceed to estimate the counting function in a passage. Here, 

, fcC^' , , 7r 2 m 2 vr 2 n 2 

a = -2> b= ^T m ' n = ~C 2 T + C 2a ik 2 ' m " ' n - 

Then, since the lattice point counting estimate is on the ellipse with semi 
axes 2 ^ and 2 ^, we have 

Nj >-ab-o{ VA + — — = A + — VA - o VA 

7T V/7T 47T 7T V 

For the upper bound we count the eigenvalues 

7r 2 m 2 7r 2 n 2 
Xm,n = + 7j2^2' m " ^ " 

Thus 

A , / /-\ 2a\/A 2b\/A kC( 1+a ^ 



A / /-\ 2a V A 26 V A A;C( 1+Q « X /fcC Q J CA /- / /-\ 
< -a&+o VA + + = A+ + — VA+o VA . 

7T \ / 7T 7T 47T V 71 " ^ / ^ ' 
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Summing, we obtain bounds for the counting function in the passages 

2M 2M r r (l+a)j C i V(P*A JH C 2i 

^ J ^ 4ir it Ait ^ tt 

j=l, even J=l> even i=Y 

v(p M ) Va c 2 -c 2 ( m+1 ) 

4vr + vr 1-C 2 ' 

where V(Pm) denotes the volume of the first M passages. We have the 
upper bound 

2M T ,,„ , 2M 



E », < ^ E 

=1, euen i=li i even x 7 

v i=i i=i / 



V(Pm) ( \ C 2a - CMl+A*) C2_ C 2(1+M) ^ 

We now are in a position to estimate the bounds of the counting function 
for the domain Q, 2 m Summing the contributions from the rooms and the 
passages, we get the lower estimate 

I n 2M L llC + C 2 - 2C 2M+1 - C 2M+2 k 2C 2a - C 2aM - c 2 °( M+1 ) \ y/X 



-A+ 



4tt \ 1-C 2 2 1-C 2a 

with the upper estimate given by 

I n 2M L hC + C 2 - 2C 2M+1 - C 2M+2 , C 2a - C2«(l+M) \ 



-A + : + £; 



a J 7T 



An \ 1-C 2 1-C 20 

As M — > oo, since there is no contribution from the tail and C < 1 we get 

j2a 

These yield the following result. 
Theorem 3.1. As A — > oo, 



VX{2C + C 2 , C 2a \ ^ AT ... \Q\ Vx {2C + C 2 W 2 



where 



Cl ^<N N (X)-W-X<C 2 ^ 
it Air it 



2C + C 2 , C 2a n „ 2C + C 2 kC 2a 
Ci = -1 — t*t ~ h — ^ > 0, C 2 = - — -=- + 



1-C 2 1 - C 2a ' 1 - C 2 1 - C 2a 

and 

C 2a 

C 2 -d = 2k —- -> as k->0. 

1 - C 2a 
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Therefore, the error is precisely of order y/X and in the limit k — > 0, we 
obtain the precise constant. 

3.2. Asymptotics for Nd(X). We now estimate bounds for the counting 
function of the Dirichlet-Laplacian. In order to do this we again look sep- 
arately at the rooms and the passages. The lower bound is obtained by 
choosing Dirichlet conditions on all the boundaries. The set of eigenvalues 
is given by 

m 2 ir 2 n 2 TT 2 /| C J ' 1 f C j for rooms, 

H — —TT- : m, n > 1 > with a = — , b = — 



4a 2 46 2 ' ' ~ J 1 ' 2 ' 2 \ kC aj for passages. 

Therefore, the contribution to the lower bound from each room or passage 
is ^A — o(\/A). Adding all contributions, we get the lower bound 

N D {\) > W.\ + o(V\). 

47T 

To get an upper bound in the first room, we choose Dirichlet conditions on 
three sides of the boundary and Neumann conditions on the right side. This 
gives us the set of eigenvalues 

(2m + l) 2 vr 2 n 2 7r 2 n \ . , t C 

+ — — : m > 0, n > 1 > with a = b= — . 



16a 2 46 2 ~ ' ~ J 2 

Counting lattice points gives ^A + + o(s/X), so as an upper estimate 
for the contribution of the first room, we get 

AC 2 C \/\ . r-, 
+ 17 +o VA). 



47T 2 7T 

For the remaining rooms and passages we choose Dirichlet conditions on the 
horizontal boundaries and Neumann conditions on the vertical ones. This 
gives us the set of eigenvalues 



m 2 7r 2 n 2 7r 2 



H — tht- '■ m > 0, n > 1 



4a 2 46 2 



Again counting lattice points gives ^A + + o(v / A). Summing all con- 
tributions yields as the upper estimate (omitting the o(\/A)-term). 





C 2M C . _ 2a l-C 2aM \ y/X 



1-C 2 2 l-C 2a J TT 

Therefore, letting M — > oo, we find the following result. 

Theorem 3.2. As X ->■ oo 

o( vi) < Wd( a) - < (|^±il + ^_ T ) ^ + 
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We remark that the lower bound is given by having Dirichlet boundary 
conditions everywhere which does not give a "v/A-term. Therefore, we can 
only get an o(vA) error estimate and not determine the sign of the \/A-term. 



4. A RELATED PROBLEM ON THE SKELETON 

We define y G dtt to be a near point of x G fHf |x — y | = dist(x, d£l), the 
distance of x to the boundary of 0. Therefore, denoting by A^(x) the set of 
near points of x, the skeleton r of is the set 

5(0) :={x£fi: card Af(x) > 1}. 

It follows that the skeleton of the R&P domain is the union of a sequence 
of line segments and parabolic arcs, T = {ej}j g N say, connecting points in 
f2; see Figure [5j For any x£ !!, there exists t G T such that x lies on one 
of two line segments C+(t),C_(t) connecting t to its 2 near points y+,y~: 
set r : O — )■ r, x — )■ t. If t(x) = t G e G T, we may therefore define the 
following co-ordinate system on r _1 (e): 

(4.1) x = x(c7, s ), r(x) = t = t(ff), a € {-1(a), 1(a)), 

where a denotes arc length along e, s is the distance from t to x along 
C+(t(a)) or C_(t(cr)), Z(<t) is the length of C + (t(a)) and C_(t(<r)) and we 
set < s < Z(t(<r)) along C+(t(<r)) and -Z(t(<r)) < s < along C_(t(<r)). 
With x = (x,y) and r(x) = (ti(x), T2(x)) it is shown in (2.4)], that the 
determinant of the Jacobian 



(4.2) J :-- 



d(x,y) 



d(a, s) 



1 



\Vr(a,s)\ ' [|Vr 1 (x)|2 + |Vr 2 (x)|2]i/2- 

For a measurable subset Tq of T and Oo := r~ 1 (ro), it then follows that, for 
any / G L 2 (Oq) with / = outside Oq; 



(4.3) / f(x.)dx = J2l da I f(<r,s) - j ds. 

Jn f-^J ej J \^r(a,s)\ 

This implies, in particular, that for / = F o r with i 7 G L 2 (Fq), 

rl(«r) 



i(<7) X 



/ 



FoWxWx = V / F(cr)(icj / — — r-ds= / F(a)a(a)da, 

o ^4 |Vr((7,s)| J Vq 



m 

where 

WW i 

(4.4) a(a) := / — — -rds. 

The integral 

(4.5) /3(a) := / \Vr(a,s)\ds 

J -1(a) 

will also feature in certain specific regions of O. 
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Figure 5. The skeleton. 



Figure 

6. r _1 (e) for 
a parabolic edge. 



We shall be considering a general R&P domain hereafter, and not the 
special case of Sections 2 and 3. Therefore we allow for the possibility that 
the Neumann Laplacian does not have a discrete spectrum. Our first task is 
to make explicit the change of co-ordinates (|4.ip in each region r _1 (e), e G T 
and then determine the map r. The edges fall into 3 groups which have to 
be handled separately. In what follows below, we denote the height of a 
room by h and of a passage by 5. 

Group 1 This consists of edges which are either in a passage or lie in the 
centre of a room with adjacent parabolic edges. Here, a = x, s = ±y, so 
the determinant of the Jacobian in (14. 2 [i equals 1 and 



1(a)- 
Hence, 

(4.6) a(a) = /3(cr 



h/2 in rooms, 
5/2 in passages. 



h in rooms, 
5 in passages. 



Group 2 These are the edges in a room which are straight line segments 
along the diagonals. In Figure [5] with the origin at O, the edge on the 
diagonal of the square given by < x < and y > 6/2 lies in this group 
and the analysis that follows is typical for all edges in this group. First, 
consider the triangle below the bisecting line, i.e. 5/2 < y < h/2 — x. 
Here, we re-parameterize points (x,y) G £1 by (cr,s), where a is the arc 
length along the skeleton measured from the corner and s is the negative 
horizontal distance of the point from the skeleton. Thus, 

(*, V)=(js + + £) , ') = (V2(£ -y),-± + x + y y 

In the triangle above the diagonal, where ^ — x < y < §,we choose s to 
be the vertical distance to the skeleton. Here, 

(x, y) = \-^=, — °-= + \ + s J , (a, s) = (V2x, -- + x + y J . 
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We note that in this whole square, we have that 

T( ^ )= (^-^)' 
with < a < Moreover, the determinant of the Jacobian J = l/y/2 

and 1(a) = a j\/2. Therefore, 
(4.7) a(a) = a, (3(a) = 2a. 

Group 3 These are edges which are such that every point on the edge 
has a re-entrant corner as one of its two near points. Thus in Figure the 
parabolic edges BE and the line segment OE belong to this group associated 
with the re-entrant corner A. We consider the parabolic edge BE. For y > 
this is determined by \AQ\ = \QQ'\, where Q = (xo,yo) is a point on the 
parabola and Q' = (xo,h/2). This gives 

5\ 2 ( b> 2 



y °-2j = { y °~2 
After a little algebra this leads to 

(4-8) y = --A- + ^ h + S ). 

In particular, the parabola intersects the x-axis 

Consider the part of the domain emanating from the re-entrant corner to 
the parabolic part of the skeleton. Let r(x,y) = (xo,yo). Then in addition 
to lying on the parabola, (xo>2/o) satisfies 

, jn . y-S/2 5 

4.9 y = —x + -. 

x 2 

The arc length along the parabola 
(4.10) c 



with to = A straightforward calculation gives that 

(4.11) a = ^{h-5) (t oy Jt% + 1 + sinh-^io) - y/2 - sinh^l)^ . 
In particular, the length of the parabolic edge is 



\CE\ = ^VWT6)-^(h-5)+^ (sinh- 1 ~ smh- 1 (1) j 

We now use (|4.8j) and (14. 9p to express to hi terms of x and y. Eliminating 
yo in (|4.8p . we get 

V- g / 2 , 6 l <u , ^ , 2 



(h-5)[ ^—^ Xo + ---(h + 8)) + x{ l = 0. 
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This yields 
(4.12) 
so that 

(4.13) 



2 2y-5 

to + — h "1 = 0, 

x 



x 



+ 1. 



For the distance s from (x,y) to (xo,yo) this then gives 



h-5 



x — 



to +[y- 



y- 5/2 h-5 



x 



2 t0 ~2 



2 _ (^ + i) 2 (x-yt ) 2 

4t§ 



and hence 



s = — 



(tg + 1) (x - Y^o) 
2t 



We need to calculate the determinant of the Jacobian ^ CT '^ . Several 

d(x,y) 

terms in this Jacobian vanish and we have 



d(a, s) 




da dto ds 


d(x,y) 




dto dy dx 



A calculation gives 

da 1 f~Z ds 

W = 2V t0 + 1{h - S) > Vx 



tj+l 
2t 



dto 
dy 



5-2y 



Therefore, 



d(a,s) 



d(x,y) 



(t 2 + lf 2 (h-5)( 



(S-2y) 



x^/(5 - 2y) 2 + Ax 2 



+ 1 



At x 



Substitution of to in terms of x and y gives 



(4.14) 



d(a, s) 



d(x,y) 



V2(h - 5) 



((5 - 2y) 2 + Ax 2 ) 



i 



(5 - 2y) 2 + Ax 2 - 5 + 2yj 



3/2 ' 



To analyse the behaviour of J 1 

5 



9{<r,s) 



near the re-entrant corner, let 



x = r cos ( 



y 



d(x,y) 

rsinO (i.e. 5 — 2y = 2r sin0). 



Then 



d(a, s) 



d(x,y) 



V2(h-5)- 



(4r 2 ) 4 



h-5 



(2r - 2rsin6>) 3/2 y/2r (1 - sin6>) 3/2 ' 

so J" 1 behaves like 1/r where r is the distance from the corner. This 
behaviour of J -1 implies that on the parabolic edge 

(4.15) a(a) < oo, but 13(a) = oo. 
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For this reason, we will need to make sure that the weight (3 does not appear 
in the analysis on those edges (like the parabolic edges) whose points have 
re-entrant near points. How we do this will be made apparent in the next 
section. 

5. A Sturm-Liouville operator 

We denote the set of edges of T which belong to groups 1 and 2 by r reg 
and those in Group 3 by r s i ng . Note that r s i ng consists of the parabolic 
edges and ones like the edge OE in Figure [5] which connect an end of a 
passage and the parabolic edges. The map r in f|4.2f) maps a re-entrant 
corner onto every point on a singular edge, which motivates us to define any 
function /one£ r s i ng to be constant. 

The underlying Hilbert spaces on T are as follows: 

(5.1) L 2 (T) = L\e;a(a)da) Cl e , 

where l e is the characteristic function of the edge e, L 2 (e; a(a)da) is the 
weighted Lebesgue space with inner-product 



(5.2) / f(a)g(a)a(a)da, 

J e 

and with / = (/ e ) 6 L 2 (T), we have f e = constant for e 6 r s i ng ; 

(5.3) H\T):= H^e-a^) Cl e , 

where for e G r reg , if 1 (e;a,/3) is the weighted Sobolev space with inner 
product 

(5.4) (f e ,9e)m(e;a,(3) = J { V/ e (a) V^)/3(<t) + f e {o)g e (a)a(a) } da. 
The inner-products on L 2 (F) and -ff 1 (r) are respectively, 

(/>ff) = E / f{°)WM<j)do 



eer ' 



and 



(/,SW) = E I {[Vf{<r)-Vg{crmcT) + [Mg{cT)]cc{a))do 

eP-T -' e 
fc-fc-L rcg 

(5.5) + l\fW)W)\^)dcJ. 

J e 



eGr sin 



We define Hr + / to be the self-adjoint operator in L 2 {T) associated 
with the H 1 ^) inner-product, where / is the identity operator on L 2 (T). 
The following theorem is readily proved by a standard argument; cf., [5], 
Theorem 3.3. 
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Theorem 5.1. The domain D(Hr) of Hp consists of sequences u = (u e ) G 
i^ 1 (r) which satisfy the following : 

(1) u e = constant for e G T s i ng ; 

(2) for each e G T reg , f3u' e is locally absolutely continuous on the interior 
of e and 

(ill 

(5.6) i im{/3((j) _^ } = o 

as a tends to the end points of e; 

(3) for e G F sing , (H r u) e = 0; 

(4) for e G T reg , 

(Hm) e (a) 1 '' 



,, ,du e 



a(a) da \_ da 

An important part in the analysis will be played by the following operator 
which maps functions on the skeleton to functions on the R&P domain. 
Define 

(5.7) T : L 2 (r) -> L 2 (fl) by T Q f = f o r for / G L 2 (T). 
Lemma 5.2. Let e G T. Then for g G L 2 (r _1 (e)) and t(a) G e we have 

(5.8) (2fa)(<7) = -L s(<7, a) - 1 da. 

i-i(«T) |Vr(<7,s)| 

Proof For F G L 2 (e), 

!{T*g)(a)F{a)da = [ g(a)(T F)(a)da 

= / Fjv) ( /"' (,T) lT f^ ,S \, da ) a((7) da, 
J e \a(a) J_ l{a) \Vr(a,s)\ J 

proving the result. □ 

Remark 5.3. Note that in particular, for t(<j) G e G r s i ng , the function 
(T^g)(a) is constant and takes the value 

(5.9) Po9)(v) = r^rrv I dx. 

\ T ( e )l Jr-i(e) 

Lemma 5.4. To is an isometry and so TqTq = /, the identity on L 2 (T). 
Proof. Let F G L(T). Then for any edge e, we have 



/ \(Fot)(x)\ 2 dx = [ \F(a)\ 2 ) ds da 

Jr-He) Je J -1(a) |Vt(<7,s)| 

(5.10) = j\F{a)\ 2 a{a) da. 

Adding the contributions from all edges shows that ||ToF|| i2 (Q) = | \F\ \ tf(r) 
which completes the proof. □ 
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The proof of the next lemma on how To interacts with derivatives is the 
same as that of (5J Lemma 3.2]. 

Lemma 5.5. For F G H 1 ^), 

(5.11) ]T I \V(For){x)\ 2 d X = £ / \F'(a)\ 2 (3(a) da. 

e&T Jr He) ee r rcg Je 

Let H\Sl) "©eer^Cr-^e)). 

Corollary 5.6. The map T% : H 1 ^) — > H 1 ^) given by T\f = for for 

f G H (r) is an isometry. 

In the final theorem, Ei denotes the natural embedding of H l {VL) into the 
space © ee r ^ 2 ( T_1 ( e ))' an d denotes the selfadjoint operator associated 
with the H l (Q) inner-product; thus Hq is the orthogonal sum 

Hq = ^(— Ajv jT - i( e )). 

Theorem 5.7. (1) The operator (H-p + I) -1 is not compact on L 2 (T) 
and G a ess (HY). 

(2) is no£ compact and G a ess (HQ). 

(3) // Ei(TiT* - I) : H 1 ^) -)• L 2 (fi) is compact, then a ess (H n ) C 
<7ess(-ffr)- 

Proof. (1) On the singular edges (#r + -0 -1 is just /. As there are 
infinitely many singular edges, (Hy + I)^ 1 is not compact. Moreover, 
by considering functions supported on one singular edge, we see that 
is an eigenvalue of Hr of infinite multiplicity. 

(2) Take an infinite sequence of edges (e n ) n£ N and consider the sequence 
of normalised characteristic functions 

^ n ( X ) = /i T7 J r-'( e „)W for x G 0. 

This bounded sequence in H 1 ^) has no convergent subsequence in 
L 2 ($7). Moreover, every ip n is an eigenfunction of Hq with eigenvalue 
0. 

(3) The proof of this is similar to that of Theorem 4.4 in [SJ. 

□ 

It is interesting to compare the results we get in Theorem 15.71 with those 
obtained in [2] for horn shaped domains where the authors are able to relate 
the essential spectrum of the Neumann Laplacian on the horn to that of a 
Schrodinger operator on the skeleton. In our case the re-entrant corners of 
the Rooms and Passages domain force us to introduce singular edges which 
cause the operator on the skeleton to decouple and prevent such a detailed 
result from being obtained. 
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